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Abstract We analyze large scale mapping observations of the molecular lines in the 
i^CO (J=2-l) , i^co (J=3-2) , i^CO (J=2-l) , and "CO (J=3-2) transition emissions 
toward the Cepheus B molecular cloud with the KOSMA 3m-telescope. The integrated 
intensity map of the ^^CO ( J=2-l) transition has shown a structure with a compact core 
and a compact ridge extended in the north-west of the core. The cloud is surrounded by an 
optically bright rim, where the radiation-driven implosion (RDI) may greatly change the 
gas properties. The intensities of the CO ( J=3-2) transition are higher than those of the 
CO ( J=2-l) transition along the rim area. We find characteristic RDI structure in position- 
velocity diagrams. Non-LTE Large velocity gradient (LVG) model analysis shows that the 
density and temperature at the edge are higher than that in the center. Our results provide 
evidences that RDI is taking place in Cepheus B molecular cloud. 
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1 INTRODUCTION 

The associati on Cepheus OB3 is one of the youngest gr oups of early ty pe stars associated with the HII 
region S155 jSharplesslI 19591) at a distance of 0.7 kpc (lBlaauwl[l964l) . The observations in the ^^CO 



( J=l-0) line have shown a few bright components named Cepheus A, B, C, D, E, F (Sargent 1977, 
1979). Of the six components, the hottest CO feature and abundant dust emission were detected in the 
Cepheus B, which was surrounded by an arc-like optical HII region S155. The molecular cloud/HII 
region complex with such a morphology is called bright rimmed cloud (BRC), where star formation 
activity can be triggered by the radiation -driven implosion (RDI) mechanism. 

The radio continuum observations jpelli et al.L Il978l) by Westerbork Synthesis Radio Telescope 
(WSRT) confirmed that the UV radiation from a bright 07 star created an arc-shaped ionization front 
which surroun ded the Cepheus B molecular cloud. The observations of the H2CO and recombina- 
tion lines (Panagia et all Il98ll) have suggested that the ionized material was flowing away at about 
11 kms^^, and the ionization front was moving into the molecular cloud at a velocity of about 2 
km s~^. ^^CO ( J=3-2) and ^^CO ( J=2-l) observations together with the far-infrared analysis showed 
that the cloud was ext ernally heated and the edge of the cloud was compressed by the expansion of S 155 
( lMinchinetaUll992h. 



Beuther et alj ( 2000h observed the Cepheus B area at a large scale in the J=3-2 and 2-1 transi- 



tions of ^^CO, ^"^CO and C^®0 with the KOSMA 3m-telescope. Based on the volume densities derived 
from FDR model, they found that Cepheus B w as highly clumped and the volume of these clumps fills 
only 2-4% of the whole cloud. iMookeriea et alj (l2006i) confirmed this with fur ther observation s of [C I] 
- '^Pi and ^^CO ( J=4-3) . Based on the studies of protoplanetary disks, iGetman I det al.h found a 
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Fig. 1 Spectra at the center positions. The vertical axis is main beam temperature Tmb, and 
the horizontal axis is the velocity. Left panel: ^^CO (J=2-l and J=3-2), right panel: ^•^CO 
( J=2-l and J=3-2). (Black for J=2-l transitions, and red for J=3-2 transitions) 



spatial-temporal gradient of young stars from the molecular core toward HD2 17086, which identified 
HD2 17086 as the primary ionizing source of the cloud. 

Although a lot of efforts have been made to characterize the physical conditions of the Cepheus B, 
the mechanism of the star formation in this region is still not quite clear We present analysis of a larger 
scale observation of multiple CO lines, and discuss the RDl as the possible trigger of the star formation 
in this region. 

2 OBSERVATIONS 

The observations were carried out with the 3m KOSMA sub-millimeter telescope using the On-The-Fly 
mode in March, 2004. Mapping observations of CO were made at ^^CO ( J=2-l) (230.538 GHz), ^^CO 
( J=3-2) (345.789 GHz), "cO ( J=2-l) (220.399 GHz), and ^^CO (J=3-2) (330.588 GHz) transitions. 
The dual-channel 230/345 GHz SIS receivers (Graf et al. 1998) were used to simultaneously observe the 
two transitions of the isotopes. The system temperatures were about 160 K for ,7 = 2 — 1 transitions and 
300K for J = 3 — 2 transitions. The pointing accuracy was better than 15" and the angular resolutions 
(FWHM of beam) were 130" at 220 GHz and 80" at 330 GHz. The efficiency of the telescope Fofr was 
93%, and the main beam efficiencies of B^s were 68% at 230GHz an d 72% at 345GH z. Line intensities 
were converted on the main beam scale, using Tmb = (-fcff/Scff)^^ (lDownesLll989l) . 

For the ^^CO ( J=2-l) and ^^CO (J=3-2) transition emissions, a 15x 15 map were observed center- 
ing at RA=22'^55"^20.708", DEC=62°20' 00.02" (B1950), while for ^^CO (J=2-l) and ^^CO (J=3-2) 
, a 12x 11-point map were made centering at RA=22'>55"'24.000% DEC=62°20' 41.01" (B1950). The 
grid spacing was 1'. Each point of the grid was observed 5-6 times with a total integration time of 2 
minutes and the spectra were averaged in order to increase the signal to noise ratio. The observed data 
were reduced and plotted with CLASS and GREG of Gildas package (Guilloteau & Lucas, 2000). The 
CO ( J=3-2) transitions were smoothed to a resolution of 130" of CO ( J=2-l) transitions. 

3 RESULTS 

The spectra of ^^CO ( J=2-l) , ^^CO ( J=3-2) , ^^CO ( J=2-l) , and ^^CO ( J=3-2) at the center positions 
of the observations are shown in Figure [T] The observed spectra show close to Gaussian line shapes. 
From the 4 lines, the average vlsr is —13.1 ± 0.1 km s^^. 

We have made integration intensity maps of the ^^CO (J=2-l) line integrated from -18 to -8 
km s^^ superimposed on the optical POSS2 red plate as shown in Figure |2] 07 star HD217086 and 
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Bl star HD2 17061 of Cephe us OB3 associat ion are marked, and the former is thought to be the primary 
ionizing source of the cloud jGetmanHet al.h . The map has shown a compact core and a compact ridge 
with a cometary shaped tail extended in south-east direction. There is a steep drop-off in CO emission at 
northern and western side of the cloud, forming a well-defined boundary edge. Meanwhile an optically 
bright rim is shown in north-west direction as seen in the POSS2 red plate in gray scale. The compact 
ridge of the CO integrated intensity contours is surrounde d by the br i ght rim . The average radial velocity 
of the Cepheus OB3 association is —14.5 ± 4.5 km s ~ ^ dOarmanv I [l973l derived from optial spectra), 
while in the north-western part the vlsr of our CO observation is also about -15 to -14 km s~^(see 
Figure|3]l. The agreement between optical and CO velocities confirms that the bright rim and the molec- 
ular cloud are physically associated. The morphology of the cloud and the HII region matches well with 
the scenario of b right-rimmed c louds, which is possibly shaped by the radiation-driven implosion of the 
lonizmg source (lBertoldiLfl989l) . 

In order to investigate the kinetic and physical conditions in this region, we take position A (offset 
1, -1) at the center of the cloud and position B (offset -2, 1) at the edge, as marked in Figure |2] The 
line properties are compared in Table [T] At position B, the ^^CO (J=3-2) transition intensity is higher 
than that of ^^CO (J=2-l) , while at position A, the ^^CO ( J=2-l) transition intensity is higher than 
that of ^^CO (J=3-2) . The reversed line ratio suggests higher density or higher temperature at position 
B. Further, we plot spectra of the ^^CO (J=3-2) (red) over ^^CO (J=2-l) (black) of all grid points in 
Figure[3] We can see that the A and B positions are not alone in this aspect. The ^■^CO ( J=3-2) intensity 
is higher than ^^CO ( J=2-l) intensity across the bright rim area for most grid positions, indicating the 
UV radiation is interacting with the molecular cloud and pumping molecular to higher excitation levels. 

Table 1 ^^CO Line Parameters 



Line 


VLSR(kms ^) 


FWHM(kms-i) 


Tn.b(K) 


^^CO (J=2-l) (A) 


-12.6 (0.1) 


2.8(0.1) 


24.2 (0.3) 


i^CO (J=3-2) (A) 


-12.9 (0.1) 


3.0 (0.1) 


21.9 (0.4) 


i^CO (J=2-l) (B) 


-13.7 (0.1) 


2.4 (0.1) 


14.8 (0.8) 


^^CO (J=3-2) (B) 


-13.6 (0.1) 


2.6 (0.1) 


19.7(1.0) 



From Figure[3]we also find a spatial shift of vlsr from the head (north-west) to the tail (south-east) 
of the cloud, which is clearly shown in position-velocity (PV) diagrams along the diagonal direction 
from north-west to south-east in Figure |4] 

4 DISCUSSION 

As seen in Figure |2] the observed molecular cloud is surrounded by an ionized bright rim facing the 
exciting source with a hot core close to the rim and a cometary shaped tail on the side opposite to 
the rim. This scenario matches well with BRCs which is modeled as externally illuminat ed, photo- 
evapo rated and ablated into elongated head-tail morphologies by UV radiation of OB stars (iReipurthi 
Il983b . The isolation and simple morphology of BRCs make them ideal environments to investigate the 
RDI (jBertoldi, 1989) mode of star formation, where pressure from the ionization shock front at the 
surface propagates through a cloud and overcomes the magnetic, turbulent and thermal pressure trigger 
cloud collapse, thereby triggering local star formation. 

In our case, the spatial shift of vlsr shows strong connections with the RDI mechanism. The PV 
diagrams in Figure |4] shows a velocity gradient from the head (north-west) of the cloud to the tail 
(south-east). The blueward elongation toward the head by about 3 km s^^from the center indicates 
that the head part is expanding towards us, while the redward elongation toward the tail by about 4 
km ^indicates the tail part is flowing away to the o pposite direction. The velo city gradients of the 
head and tail matches with the RDI model prediction of lLefloch & LazareffI (1 19941) . where the red-shift 
of the tail is caused by the acceleration of the gas across the tail, and the blue-shift of the head is because 
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Fig. 2 Contours of the ^^CO ( J=2-l) line integrated intensity from -18 to -8 km s~^(white 
solid line) overlaid on POSS2 Red image (gray scale). The contour levels of solid lines start 
from 11 to 70 by step of 7.4 Kkm s^^. Points A and B are marked as typical positions at the 
center (A) and the edge (B) of the cloud for further analysis. The bright rim in optical image 
is marked in yellow. HD217086 and HD217061 are marked in orange cross. 



the compressed gas in the head moves faster than the gas immediately behind it. Such PV diagrams 
are reported by various BRC and comet ary globule studies dLefloch & Lazare5ll995l : ISugitaniL 19971 : 
I White et al.Lfl997l:lBachiller et al.Ll2002h which are explained by the RDI mechanism, and so are taken 
as a characteristic structure of the RDI. 

The reversed line ratios across the bright rim are another interesting feature. The spatial distribution 
of the reversed line ratios suggest that they are physically connected with the RDI along the rim area. 
Higher transitions are more readily excited by higher temperatures or higher densities (Petitpas, 1998), 
either of which can be a result of RDI. The highly luminous ionizing shock front from the RDI interacts 
with the molecular cloud, pumping more particles from a lower level to a higher one, which makes the 
J=3-2 emission stronger along the rim area. 
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Fig. 3 Spectra of the ^^CO ( J=2-l) (black) and ^^CO ( J=3-2) (red) of all grid points. For 
each grid point, the horizontal axis is the velocity with a scale from -30 to 5km s~^, and the 
vertical axis is the main beam temperature from -2 to 25 K. A and B are the same positions 
as marked in Figure |2] 
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Fig. 4 P-V diagram constructed from the ^^CO ( J=2-l) (left panel) and ^^CO ( J=3-2) (right 
panel) along the diagonal direction from north-west (positive offsets) to south-east (negative 
offsets). 



In order to investigate RDI effects and i ts physical characteristics quantitatively , we made further 
analysis with a large velocity gradient (LVG. lGoldreich et al.lll974l : IOin et al l l2008ah radiative transfer 
model in the MIRJAD software package. In this model, the molecular cloud is considered spherically 
uniform and the non-LTE excitation properties are calculated based on the kinetic temperature (Tkin), 
the column density per unit velocity interval (Nco/AV), and the H2 density (nHa)- In this case, we use 
the line intensity of ^^CO (J=2-l) and the hne ratio of ^^CO (J=3-2) /^^CO (J=2-l) to determine the 
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Fig. 5 LVG model analysis at position A (left panel, Nco = 4.6 x lO^^cm"^) and B (right 
panel, Nco = 1-2 x W^cmr'^) using the ^^CO (J=2-l) and ^^CO ( J=3-2) line. The hor- 
izontal and vertical axes stand for volume density n(H2) and kinetic temperature Tkin in log 
scale, respectively. Red curves are line intensity of ^^CO ( J=2-l) and blue curves are the line 
ratios of ^^CO ( J=3-2) /^^CO (J=2-l) . The results of Tkin and n(H2) are taken from the 
cross points of certain red and blue curves. The error bars indicate the range of the possible 
solutions listed in Table |2] 

volume density nna and kinetic temperature T^in at the two selected positions A and B, assuming the 
column density in a range of 1 x 10^"' to 1 x lO^^cm^^, and AV=1 km s^^. Representative solutions 
that can be fitted to the observations are listed in Table|2] A typical case is shown in Figure[5] 

Table 2 Solutions of LVG analysis 



Position 


Nco (10^'' cm-2) 


n(H2) (10* 


Tkin (K) 


A 


3.2 


1.6 


40 




6.8 


3.2 


31 




10 


4.0 


31 




22 


2.0 


30 




46 


1.0 


30 




68 


0.63 


30 




100 


0.20 


32 


B 


1.1 


4.0 


71 




1.2 


8.9 


50 




1.3 


50 


45 



For position A, various column densities are possible to explain the observation results self- 
consistently covering a range from 3 x 10^® to 1 x VS^^cjmT^, with corresponding volume densities 
ranging from 0.2 x 10"' to 4.0 x lO'^cm"'^. The average value of the possible volume densities is 
2.0 X lO^cm^'^. For position B, the few self-consistent solutions require the column densities in a 
tiny range of 1.1 x 10^^ to 1.3 x VS^^cmT^ . The average density derived is 1.5 x Vd^cmT^, nearly one 
order of magnitude higher than that at position A. The average kinematic temperature of B (61 K) is 
also twice as high as that of A (32 K). Even in the cases that both positions have similar densities (e.g. 
n{H2) = 4.0 X 10'' cm^'^), the kinematic temperature required at position B is still 2.3 times than that 
at position A. Therefore we can conclude that the reversed line ratios along the bright rim are produced 
by either high densities, or high temperatures, and likely both. Both possibilities can be explained by 
radiation and pressure of the ongoing RDl, where ionizing shock wave of the RDl propagates into the 
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molecular cloud, compressing it at the edge, and thus makes the density and temperature of the edge 
even higher than the center. 

In comparison with the LVG analysis results, we have also calculate d the column densi ties assuming 
LTE conditions. Taking th e upper energy leve l temperature of 16.6 K jOin et al l l2008bh . the column 
density of CO is given by dGarden etaUll99l[) : 



and 



where the beam filling factor / is assumed to be 1, and the cosmic background temperature T^g is taken 
as 2.732 K. With the ^^CO ( J=2-l) line properties listed in Table[T] we derive a CO column density of 
4.0 X 10^^cm~^ for position A, and 1.9 x lO^^cm"^ for position B. Because ^^CO (J=2-l) is always 
considered optically thick, the column density derived here should be a lower limit. Compared with the 
column densities from LVG analysis, we find that LVG results at position A roughly agree with the LTE 
calculation, but for position B, all the LVG solutions are lower than the LTE lower limit. This result 
agrees with the scenario that at position B the LVG calculation requires less column density than LTE 
to explain the observing results because it takes into account enhanced photon escape probability due 
to greater velocity gradient at the bright rim, where the nearby HII region agitates the molecular cloud 
with strong FUV radiation pressure and heating, exhibiting a case of RDI. 



5 SUMMARY 

We present analysis of large scale mapping observations of the molecular lines ^^CO (J=2-l) ^^CO 
(J=3-2) ^■^CO (J=2-l) and ^'^CO ( J=3-2) transition emissions toward the Cepheus B molecular cloud. 
The target lies at the edge of HII region S155, and is surrounded by an optically bright rim. The CO 
J=3-2 emission is stronger than J=2-l emission along the bright rim. PV diagrams of this region show 
a velocity gradient from the head of the cloud to the tail, which is characteristic of RDI. According to 
the LVG analysis, the density and temperature at the edge of the cloud are higher than that at the center. 
All of our results provide evidences that RDI is taking place in Cepheus B molecular cloud. 
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